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1. Overview 

Using a multi-plane lensing method that we have developed (Premadi 1996; Premadi, Martel, & Matzner 
1998), we follow the evolution of light beams as they propagate through inhomogeneous universes. We 
use a P 3 M code to simulate the formation and evolution of large-scale structure. The resolution of the 
simulations is increased to sub-Megaparsec scales by using a Monte Carlo method to locate galaxies inside 
the computational volume according to the underlying particle distribution. The galaxies are approximated 
by isothermal spheres, with each morphological type having its own distribution of masses and core radii 
(Jaroszynski 1992). The morphological types are chosen in order to reproduce the observed morphology- 
density relation (Dressier 1980). This algorithm has an effective resolution of 9 orders of magnitudes in 
length, from the size of superclusters down to the core radii of the smallest galaxies. 

We consider cold dark matter models normalized to COBE, and perform a large parameter survey by 
varying the cosmological parameters S1 j Ao, H , and n (the tilt of the primordial power spectrum). The 
values of n are chosen by imposing particular values or og, the rms mass fluctuation at a scale of 8/i _1 Mpc. 
We use the power spectrum given in Bunn & White (1997). Table 1 gives the values of the parameters for 
all models. This is the largest parameter survey ever done is this field. 
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Table 1. Parameters for the various models 
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2. The Ray Shooting Experiments 



For each model, we perform numerous ray-tracing experiments, propagating beams of 341 2 light rays 
back in time up to redshifts z = 3 or 5. Distance between two neighboring rays is about 1 arcsecond. This 
work is still in progress. So far, we have performed between 20 and 40 experiments for each model. Even 
with these small numbers of runs, we can already see significant differences between the various cosmological 
models. Figure 1 shows the individual contribution of each lens plane to the shear, as a function of the 
lens redshift, for sources located at z = 3 and z = 5, and for 3 cosmological models: An Einstein-de Sitter 
model, an open model with flo = 0.2 and Ao = 0, and a flat model with flo = 0.2 and Ao = 0.8 (these results 
are obtained by averaging over many experiments). The lenses planes contributing the most are located at 
intermediate redshifts, in spite of the fact that structures are more evolved at low redshifts, and galaxies are 
more crowded at high redshift. The shear peaks at lower redshift for the Einstein-de Sitter model than for the 
other two models. We interpret this effect as a consequence of the freeze-out of the growth of perturbations 
in models with f2 < 1. In such models, structures present at low redshift must have formed earlier than in 
the Einstein-de Sitter model. Hence, such models tend to have more structure at intermediate redshift for a 
given value of as- 




Fig 1. Individual contribution of each lens plane to the shear vs redshift of the lens plane, for sources located 
at z — 3 (left panel) and z = 5 (right model). Solid curves: Einstein-de Sitter model; dotted curves: open 
model; dashed curves: flat, cosmological constant model. 

Using these experiments, we can locate actual sources on the source plane, and compute their images 
on the image plane. In Figure 2, we plot the images of distant circular sources, to illustrate interesting cases 
that occurred in some experiments: magnification and shear (Fig. 2a), double image (Fig. 2b), Einstein ring 
(Fig. 2c), and triple image (Fig. 2d). 



3. Preliminary Results 

We can use the properties of the images to compute lensing statistics that can eventually be compared 
with observations. This is the ultimate goal of this project, and it will require many more experiments 
than the ones we have performed so far. For illustrative purposes, we computed, for various models, the 
probability that a lensed quasar will have multiple images. The results arc presented in Table 2. As we see, 
the probabilities are sm all. Hence, most lensed quas ar s arc simply magnified, with no image splitting. The 




Fig. 2. Images of distant circular sources 



Table 2. Multiple- image Probabilities 
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